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On the basis of Bethe ansatz solution of one dimensional Kondo model with electronic interaction, 
the thermodynamics equilibrium of the system in finite temperature is studied in terms of the 
strategy of Yang and Yang. The string hypothesis in the spin rapidity is discussed extensively. The 
thermodynamics quantities, such as specific heat and magnetic susceptibility, are obtained. 
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I. INTRODUCTION 

It is known that the study of exact solutions is help- 
ful for the understanding of non-perturbative effects in 
the strongly correlated electronic systems. The exact so- 
lution of the one dimensional Kondo [Q model with lin- 
earized dispersion in the absence of electronic interaction 
was found in (^^] . The model with quadratic dispersion 
in the presence of electronic was shown to be exactly 
solvable at some value of electron-impurity coupling Q . 

In present paper, we study the thermodynamics of 
Kondo model in the presence of the electronic interac- 
tions. The general thermal equilibrium are discussed ex- 
actly on the basis of the known Bethe ansatz solutions 
of the model 0]. The specific heat and magnetic sus- 
ceptibility are obtained analytically in general and given 
explicitly in strong-coupling limit. The specific heat and 
the magnetic susceptibility at low temperature are dis- 
cussed. Next section we briefly exhibit the model under 
consideration and its Bethe ansatz solution. In section 
III, we demonstrate the string hypothesis and write out 
the Bethe ansatz equation in the presence of complex 
roots. In section IV, we consider the thermodynamics 
limit by introducing the density of roots and holes. In 
section V we derive the free energy of the system at ther- 
mal equilibrium according to the strategy of Yang and 
Yang ||. In section VI, We calculate the thermodynam- 
ics quantities, such as specific heat and magnetic suscep- 
tibility. In the case of strong coupling limit, we find that 
the contributions of both electrons and impurity to the 



specific heat and the magnetic susceptibility is Fermi- 
liquid like. 



II. THE MODEL AND ITS SPECTRUM 

The model Hamiltonian of a correlated electronic sys- 
tem which we consider reads 

H Q = J2s(k)C£ a C ka 

k 

+ + k2 ' k 3 + k i) C LaGLa C k 2 aC kia , 

ki,k 2 ,k3,k 4 

where Cka annihilates an electron with momentum k and 
spin component a, and e(k) = k 2 /2 (in units of % and of 
the electron mass). The electrons are coupled by both 
spin and charge interactions to a localized impurity, 

Hi = j*:(o)s Qb * 6 (o) • s° + y*;(o)* a (o), 

where the field ^> a is the Fourier transform of Cka, S° is 
the spin of the impurity and S is the spin of the electrons 
in the band. 

The present model was solved exactly by Bethe ansatz 
with periodic boundary conditions H) > It was shown that 
the model is integrable when the correlation strength u is 
proportional to the J, the strength of electron-impurity 
coupling via spin. The Bethe Ansatz equations for the 
spectrum are 
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where 9(kj 



2 tan 1 (kj/u). In the approximation 



kj ~ fc; for any j, I, the S matrix of electron-electron 
will be independent of u. Then the Yang-Baxter equa- 
tion will give no restrictions on between u and J. This 
makes it easy to understand the usual Kondo problem 
where the linear dispersion relation is adopted, whence 
the model is solvable at any value of J. 



III. STRING HYPOTHESIS 



the Bethe ansatz equation ||. The complex roots with 
the same real part A^ form a n-string, 



Ap m = Xp + i-m + 0(exp(— <WV)), (S > 0) 
m = —n + 1, — n + 3, • • • , n — 3, n — 1. 



(2) 



The set of roots {X v \v = 1,2, - ■ ■ , M} is then partitioned 
into a set of n-strings {Ao|m = — n + 1, — n + 3, • • • , n — 
3, n - 1; ft = 1, 2, • • • , M n }. Obviously, 



For the ground state (i.e. at zero temperature), the 
k's and A's are real roots of the Bethe ansatz equation 
(|l|). For the excitated state (i.e. at nonzero tempera- 
ture), however, they can be complex roots HQ. We will 
not take account of the complex roots in the charge sec- 
tor k for repulsive interaction since it does not happen 
at low temperature. The complex roots A in spin sec- 
tor always for a "bound state" with several other A's, 
which arises from the consistency of both hand side of 



M 



^ nM n , 



71=1 



where M n denotes the number of n-strings. 

Substituting those n-strings into the eq. ([!]) , we can find 
that the product of the fractions for the roots within the 
same n-string reduce to (X^—kj+inu/2)/ (X^—kj— inu/2) 
because of the alternative elimination between denomi- 
nator and numerator. Hence the Bethe ansatz equation 
(TO) becomes 



\ n — k 

-ikjL _ g-ifHfcj) J J '_ 



- inu/2 



XV, — kn — inu/2 ' 

n P 3 1 



and 



^ - ^ ^ t u _ Agp _ iu /2 * Aff> -h- iu/2 
11 11 A^-AS'-iu A7 + iu/2j-\A7-k l+i u/2- 

j3m q= — m+l P ' 1=1 ' 

The product of those equations (|J) for p = —n + 1, — n + 3, • • • , n — 3, n — 1 gives rise to 

- K+i(m + n)u/2 

( i) n TT x 

v ; 11 X m_ X n_ i(rn + n)u/2 



fim 



A™ - A™ + i(m + n - 2)u/2 Xf-Xl + i(\m - n\ + 2)u/2 
A™ - AS - i(m + n - 2)u/2 " ' XJJ - X n a - i(\m - n\ + 2)u/2 

~ K + i\™ - n\u/2 _ Ag - inu/2 ^ Ag - h - inu/2 

inu/2 

j i ^2®n/2{Xp -kj), 



: TT - 

A ™ ~ K ~ A m - "I"/ 2 K + ™/2 Y A2 - h 
Taking the logarithm of eq.Q and eq.(||) we have 

2tt _ 1 
L 



fin 
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(->„ 2( a;: ) - V e n/2 (x n a - h) = 2^ - ^ A„ mp e p/2 (A^ - as), 

f3mp 



where @ p {x) — 2 tan and 



.4 
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1, for p = m + n, \m — n|(^ 0), 

2, for p = n + m — 2, n + m — 4, • • • , |n — n| + 2, 
0, otherwise. 



(3) 



(4) 



(5) 
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The Ij and J™ are quantum numbers, the Ij are integers 
or half-integers depending on whether M is even or odd, 
the J™ are integers or half-integers depending on whether 
N — M„ — n + 1 is even or odd. 



IV. THE THERMODYNAMICS LIMIT 

The trasendental equations @ for the real parts of the 
complex roots are obviously difficult to solve. It will be 



convenient to consider the thermodynamics limit ||, i.e. 
N — > oo, L — > oo but D = N/L is fixed. Introducing the 
density distribution of roots and holes 



(7) 



we obtain from (g) the following set of integral equations, 



p(k) + p\k) = i- - i^(fc) + £ ^„/ 2 (fc|A')a(A'), 

n=l 

-„(A) + a£(A) = y^„/ 2 (A) + i^„ /2 (A|A : , Kfc , )-E^n™ ^) ^ P /2(A|A'K n (A , ), 



(8) 



nip 



where K n (x) = it 1 nu/(n 2 u 2 + x 2 ). We have adopted 
a notation convention K n (x\y)p{y) = / A" n (a; — y)p(y)dy 
etc. in the above. 

In term of the density distributions, the energy and 
the concentration of electrons as well as the number of 
down spins are given by 

/oo 
dkp(k)k 2 , 
-OO 

n r°° 

D = - = J dkp{k), 
-y = ^2 71 \ d\<j n (\). 

L n=l 

Thus the magnetic moment of the system is, 

1 

2 
L 
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M = -{N-2M) + Sl mp (9) 
p(k)dk -L^n / a„(A)dA + ^ mp , (10) 

-OO i J —QO 



n=l 

where S? mp stands for the spin of the impurity, and the 
g factor is put to unit. 

The ground state of the present model is a Fermi sea 
described by p{k) with real rapidity A, i.e. p(k) = for 
\k\ > k F and p h {k) = for |fc| < k F ; oi(A) ^ but 

(A) = (n > 1), which is the case at zero tempera- 
ture. Away from zero temperature, The density distribu- 
tions of roots and holes with respective to the momentum 
k and the spin rapidity A should be determined by the 



principles of statistical physics. Next section we will dis- 
cuss this issue extensively on the basis of the strategy of 
Yang and Yang || . 



V. THERMAL EQUILIBRIUM 

For a given p(k) and p h {k), the number of roots and 
that of holes in the neighborhood dk are Lpdk and 
Lp h dk respectively. Obviously, the total number of roots 
and holes in the neighborhood is L(p + p h )dk. For a 
given <7„(A) and tr^(A), La n dX and La^dX give rise to 
the number of n-string and the number of the vacan- 
cies of n-strings (holes) in the neighborhood d\, while 
L(a n + <Jn)d\ gives rise to the total number of n-string 
and vacancies of rt-string. Thus the total number of the 
possible choice of state in dkdX being consistent with 
given distribution functions in both charge and spin sec- 
tors is 



H(fc,A) 



[L(p + p h )dk]\ -p,- [L(a n + a*)d\]\ 



n 



[Lpdk]\[Lp h dk] \ 11 [La n dX\\[La%d\] 



As the total number of all possible state for given distri- 
bution functions is 



5 = I]s(M), 



the total entropy S will be obtained by taking logarithm 
of S, 
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S/L = J {[p(k) + p h {k)} ln[p(k) + p h {k)} - p(k) lnp(fc) - p h (k) \np h (k)} dk 
I { K (A) + o£(A)] lnK(A) + ttf (A)] - a n ( A) In a n ( A) - a h n ( A) In a h n ( A) } dX, 



(11) 



where the Boltzmann constant is put to unit. 

In the presence of the external magnetic field H, we 
must add a Zeeman term to the original Hamiltonian. 
As the Zeeman term commute with the original Hamil- 
tonian, the Bethe ansatz solution is still valid for present 
case. Therefore the energy of the system in the presence 
of external magnetic field will be 

E/L = I (k 2 - H)p{k)dk + j>2 2nH f <r n (k)dk, (12) 

J 71=1 J 

In order to obtain the thermal equilibrium at tempera- 
ture T, we should maximize the the contribution to par- 
tition function from the state described by the density 



distributions of roots and holes. As maximizing the par- 
tition function is equivalent to minimizing the free energy 
F = E-TS-pN. Here S and E are given by eq.(|ll|) and 
eq.(|l^), p is the chemical potential for canonical ensem- 
bles. The p plays the role of the Lagrangian multiplier 
for the condition L J p(k)dk = N = constant if one min- 
imizing the Helmholtz free energy = E — TS. This 
constranit implies that the essemble has fixed number of 
partiles. A constraint that the number of down spins are 
fixed was imposed in ref. || when discussing delta Fermi 
gas, whereas we will not impose no physics constraints in 
the following discussion. 

Making use of the relations derived from eq. (0) 



5p h (k) = -5p(k)+Y,K n/2 (k\\)6<j n (\), 

n 

5a h n {X) = -5a n {\) + K n/2 (X\k)6p(k) ~ ^ A mnp K p/2 (X\X')Sa m (X'), (13) 

mp 

we obtain the following conditions from the minimum condition §F — 0, namely 

e(k) = -p + k 2 -H- T^K n/2 (k\\) ln(l + e -^' T ), (14) 

n 

f„(A) = 2nH - TK n/2 (X\k) ln(l + er< k ^ T ) 

+TY / A nmp K p/2 (X\X') ln(l + e-«™( y )/ T ), (15) 
mp 



where we have written 
P h {k) 
p{k) 

°"n(A) _ 

o-„(A) 



exp[e(fc)/T], 



exp[£ n (A)/T]. 



Principally, once e(k) and £(A) are solved from eq.([l4j- 
|l5| ), the equilibrium distributions p(k) and a n {X) at tem- 
perature T will be known from the following relations, 



p(k)[l + e^/ T } = 



1 1 

2^ ~ I 



K 1 (k) + ^2K n/2 (k\X')a n (X'), 



a n (X)[l + e« A )/ T ] = \K nl2 {X) + K n/2 {X\k')p{k') 

— ^ A nmq K q / 2 (X\X')a m {X'). 



(16) 
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(17) 



The free energy per unit length reads 

F/L = J dkp{k) [fc 2 - e{k) -H-T(l + e £(fe)/T ) ln(l + e - £ «/ T 
+ J2 I d\a n (X) [2ni/-6 l (A)-T(l + e«" (A)/T )ln(l + e-«" (A)/T ) 

Integrating eq.(|l4j) over A; after multiplying it with D~ 1 p, we get the chemical potential 
H = i J (fc 2 - e(k) - H)p{k)dk 

H / / K ^ k - A ) ln (! + e-^ (x)/T )p(k)d\dk. (18) 

n J 

Integrating eq.(|l5|) over A and summing over n after multiplying it with D~ 1 a ni we obtain the following relation 

J2 ! in{X)cx n {\)dX = ^2nH I a n (X)dX 
-TJ2[ [ K n/2 (X-k)ln(l + e-<W T )a n (X)dkdX 

n J 

+T A ram I I K q/2 (X - A') ln(l + ( A ' )/T )a n ( X)dX'dX. (19) 



Using the relations (|18[jl9|) , we can write out the free energy in terms of e and £ only, 



F = [iN + Tj (K 1 (k)-^-)ln(l + , 



-e(k)/T 



)dk 



-T W K n/2 (X) ln(l + e-^/ T )dX. (20) 

Consequently, the partition function is given by 

Z = e- F / T , (21) 

where the Boltzmann constant is put to unit. 

The thermodynamics functions, partition function Z, 
free energy F, and thermal potential f2 etc, are of impor- 
tance, as knowing either of them, one is able to calculate 
all thermodynamics properties of the system in principle. 
However, it is difficult to obtain an analytic expressions 
of e(k) and £ n (A) from the coupled non- linear integral 
equations (|l4| ) and ( |l5| ) . So we are not able to derive ex- 
plicit results for thermodynamics quantities. Moreover, 
we will obtain some plausible results for some special 
cases in next section. 



VI. THERMODYNAMIC QUANTITIES 

In general, the free energy of our model should be cal- 
culated by formula (p0| ), where the e(k) and £(A) are de- 
termined by eq.( ^4[]l5| ). Then the other thermodynamic 
quantities are obtainable from thermodynamic relations. 



In the Appendix, we show that if p, e and are 
implicit functions of some thermodynamic quantities x 
(such as T, L), the derivative of eq.(|2(]) with respect to 
x is the same as the partial derivative of eq. (j2^) with 
respect to the explicit variable x. It is easy to get the 
pressure of the system 

P = -JL = ^ /Ml + e- eW/T )*, (22) 

which is the same as Yang and Yang's expression for- 
mally. However, the integral equation which e(k) obeys 
are different. In our present case the contributions from 
both impurity and the spin rapidity are involved. 

In terms of e and £, the entropy S — —(dF/dT) be- 
comes 



S= - J ln(l 



+ e -m/T 
e -m/T e(fc) 

l + e -e(k)/T T 



L 

2^ " 
L 

~27T 



1 



7T U 
1 



-fc 2 

u 



dk 



nu/2 



iru 2 + k 2 
ln(l + e-^/ T )dX 



dk 



E 



7T (nu/2) 2 + A 2 

1 nu/2 £ n (\)/T 
Tr(nu/2) 2 + X 2 l + e^W/ T 



dX. 



(23) 



The other thermodynamics quantities is formally obtain- 
able, e.g. 
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For the sake of saving space, we did not write them out. 

The free energy (^2|) is conveniently partitioned as two 
parts. F — F e + Fi, here Fi stands for the contribution 
from the impurity, 



± i — 



7T J U 2 + k 2 

nu/2 

7T ^ J (nu/2) 2 ri 2 + A 2 
and F e for that from the electrons 



ln(l + e- e ^/ T )dk 



A. Strong coupling limit 

For the strong coupling u >> 1 and non- vanishing ex- 
ternal magnetic field, 

we keep the leading term in cq.([l4|) and (|l5|), 



e(k) — k 2 - H - fi, 
£»(A) = 2nif. 



(25) 



r, =r / ^^i^in[i + e( f '- fe2 )/ T ] 



7T U 2 + fc 2 

dA mt/2 



7T (nw/2) 2 + A 2 



ln[l 



,-2nH/T ] 



L 



F e — \iN — T— / ln[l + e("- fe ^ T }dk, (26) 

27T J 

where rj = /i + H . Since tan _1 (fc/u) ~ k/u for u » 1, 
integrating eq.(p6|) we have 



Fa 



ttw 7 1 + e ( fe2 -")/ T 



2L f°° k 2 
F e = fiN / — — —dk 



dk-TY,H± + e~ 2nH/T ), 
(27) 



7T J Q l + e ( k2 -v)/T 
Now we compute the common integration in eq. (|2 



1=1 7 7777- T—T^Tcdk, 



l Q 1 + e (k^eta/T) """> 

Changing variables of z = (k 2 — f])/T brings it to 
2 J-n/T 1 + e z 



(28) 



(29) 



The free energy related to impurity and electrons become which is conveniently slit up into three terms, 



T r i/T , T r l/T (n - 



T r l/T (n-zTfl 2 , T ' v 
dz + — 



(rj + zT) 1 / 2 



1 + e 2 



dz. 



(30) 



The right end of the interval for integration in the second term can be regarded as oo due to the contribution from 
large value of z is negligible. Integrating eq. (pfl) after expanding the numerator as Taylor series, we obtain 



71 = 



n+l/2 



1+1/2 



(31) 



where £(2n + 2) is the Ricmann Zeta function. Conse- 
quently, the free energy is obtained 



Fi=— I-TY\n[l + e- 2nH ' T }, 

IT 11 ' 



2L 

F e = fiN 1. 

IT 



(32) 



B. Thermodynamic quantities at low temperature 

In the low temperature approximation, eq.(|3^) be- 
comes 



-i-ry, 



-2nH/T 



(33) 



The magnetization: 

The contributions of electrons and the impurity to the 
magnetization are obtained 



2L 

M e = —I h , 
7T 

4 

nu 



(34) 



n=l 



where 



dl d 2 I d 2 I 

h = m , ^ = ^2> ht = W (35) 



G 



Clearly, the system has spontaneous magnetization, 
L i in ttL . 



M e {H^0) = -n 



1/2 



24 



tV 3/2 , 



7TU 12li 

Specific heat and magnetic susceptibility: 
Specific heat and the magnetic susceptibility are 



Ce - 


2LT 

Itt, 

TT 








Xe = 


Ihh, 

TT 








Ci = 


4T 


AH 2 


e -2H/T 




Itt + 

TTU 






Xi = 


4 

hh 4 

TTU 


4 

-T e ~ 


-2H/T 


(36) 



where C e {Ci) is the specific heat of electrons (impurity), 
and Xe(Xi) i s the magnetic susceptibility of electrons (im- 
purity). 

Wilson's treatment of Kondo model by re- 
normalization-group calculation has made it possible 
for the determination of proportionality factor( "Wilson 
Ratio") relating low-temperature and high-temperature 
dimensional scales. When T — > 0, because e~ 2H / T de- 
crease more rapidly than T, the second term of eq.(|36|) 
can be neglected. Then we are able to evaluate the 
Wilson's Ratio 



R = 



Xi/Xe 
O/CTe 



1. 



Furthermore, if we ignore the small term in I u and Ihh, 
that is we only remain the first term of them, we find 
that the impurity's contribution to specific heat at low 
temperature is Fermi-liquid like 



H) -l/2 T 



(37) 



and so is the magnetic susceptibility 

irT 2 



Xi 



■-!-&.■ 

TTU 



H 



-1/2 



8u 



(/i + Jf)- 6 / 2 . (38) 



Obviously, the zero temperature susceptibility is of finite 

1 



x*(r = o) 



-Qi + H) 



'1/2 



(39) 



indicating that the impurity spin manifest in high tem- 
perature regime by Curie's Law, Xi 1/T, is now com- 
pletely screened. We interpret the effect as due to strong 



effective coupling between impurity and the conduction 
electrons leading to the formation of a local singlet, and 
the infra-red physics is dominated by a strong coupling 
fixed point. 

The electron's contributions to specific heat and mag- 
netic susceptibility are 

Xe = ^0* + ^ + + H)-^. (40) 

We have analyzed the thermodynamics of Kondo 
Model with electronic interactions, particularly discussed 
the case of strong coupling limit extensively. In that case 
we have shown the impurity's contribution to the specific 
heat and magnetic susceptibility of the system is Fermi- 
liquid like and shown that in very low temperature the 
system has the property of spontaneous magnetization. 



APPENDIX 



Since e and £ which solve eq( |14[]15| ) evidently depend 
on (j,, we should consider them as functionals of \i which 
is usually functions of some thermodynamic variables x 
(such as T or L). The derivative of free energy (20) is 
given by 



OF 

dx 



dk 



N 



L/2-k- 



dft 

K x {k) de dfi 



1 + eA k )/ T dpi dx 



1 



•W/ T dfi dx' 



(41) 



where (dF/dx)^^^ stands for partial derivative of F with 
respect to the explicit variable x while fi, e, and £ are re- 
garded as irrelevant to x. 

The derivative of (|l4|) with respect to pt is 



de 
9/i 



E 



dX 



K n /2(k - A) d^n 



(42) 



Integrating ([ 
p, we obtain 



1 + e €n(A)/T dfl 
ft after multiplying both hand sides with 



de 

dk-p(k) 



4 + E/^(A)-iE/^ T 



E 

nmq 1 



d\d\ 



i d£, n A Tl 



(43) 
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In deriving to the above equation, the second equation of eq. (|16|) has been used. We take derivative of eq.(|15|), then 
multiply both hand sides with <j n and integrate over A. Summing over the subscript n in what we obtained, we have 

9 £n ? \ \ f f ji n K n/2(^ - k)a n (X) 



/ d\^*n(\)=Yl J J dMx 



1 + e e ( fe )/ T 



D / / " A '| i + ^)/t V(A - A' K (A). (44) 



nmg 

With the help of the first equation of eq. (|l6|), eq.([44"|) and eq.(|43|) give rise to 

Z/27T ^ y # n / 2 (A) 



N+ dk 



l + e e(k)/T dfl 



TfdX T < J /T ^ = 0- (45) 



Thus the complete cancellation of the last three terms in 
eq. ( fTl| ) concludes that 



dF 
dx 



dF 

dx 
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